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Abstract 

Starting from the Rodrigues representation of polynomial solutions 
of the general hypergeometric-type differential equation complemen- 
tary polynomials are constructed using a natural method. Among the 
key results is a generating function in closed form leading to short 
and transparent derivations of recursion relations and addition theo- 
rem. The complementary polynomials satisfy a hypergeometric-type 
differential equation themselves, have a three-term recursion among 
others and obey Rodrigues formulas. Applications to the classical 
polynomials are given. 
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1 Introduction 



Real polynomial solutions Pi(x) of the hypergeometric-type differential equa- 
tion (ODE) 

(P P riP J 

a{x) id + T{x) ii + AlPl{x) = °' A/ = ~ y ~ 2 (/ " l)a " (1) 

with 1 = 0,1,... and real, first and second-order coefficient polynomials 

a(x) = ex 2 + 2fx + g, r = a\ + b\x (2) 

are analyzed in ref. [Tj, [2J. The (unnormalized) polynomials are generated 
from the Rodrigues formula 

Pi{x) = -^^-Aa l (x)w(x)), 1 = 0,1,..., (3) 
w(x) dx l 

where w(x) is the possibly I dependent weight function on the fundamental 
interval (a, b) that satisfies Pearson's ODE 

a(x)w'(x) = [t(x) — a'(x)]w(x) (4) 

to assure the self-adjointness of the differential operator of the hypergeomet- 
ric ODE. Polynomial solutions of ODEs with I dependent coefficients are 
studied in ref. [3] along with their orthogonality properties and zero distri- 
butions, which we therefore do not address here. 

Here our first goal is to construct complementary polynomials for them 
by reworking their Rodrigues representation, Eq. ([3]), in a simple and nat- 
ural way. The generating function of these complementary polynomials is 
obtained in closed form allowing for short and transparent derivations of 
general properties shared by the complementary polynomials. 

The paper is organized as follows. In the next section we introduce and 
construct the complementary polynomials. In Section 3 we establish their 
generating function, the key result from which recursion relations and an 
addition theorem are derived in Section 4. The Sturm-Liouville ODE is 
derived in Section 5. Classical polynomial examples are given in Section 6. 
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2 Complementary Polynomials 

Definition. We now introduce the complementary polynomials V u (x;l) 
defining them in terms of the generalized Rodrigues representation 

PA*) = (<r{xT- v w{*yP»{x\ 0) , (5) 



where v = 0, 1, . . . , I; I = 0,1, ... . 

Theorem 1. V u (x; I) is a polynomial of degree v that satisfies the recur- 
sive differential equation: 

dV (t- 1) 

V u+1 (x; I) = a(x) ^' ] + [r{x) + (I — u — l)a'(x)]V v (x; I). (6) 

By the Rodrigues formula (j3J), Vq{x;1) = 1. 

Proof. Equations ([S]), and ([SD follow by induction. The first step, v = 1, 
is derived by carrying out explicitly the innermost differentiation in Eq. ([3]), 
which is a natural way of working with the Rodrigues formula ([3]) that yields 

p^x) = ] f / (la l - 1 (x)w(x)a'(x) + a l (x)w'(x)) (7) 
w(x) aar 1 v ' 

showing, upon substituting Pearson's ODE (j4j), that 

Vi(x;l) = (l-l)a'(x)+r(x). (8) 

Assuming the validity of the Rodrigues formula (jSJ) for v we carry out another 
differentiation in Eq. (J2D obtaining 

1 ( 



P ^ = dx i-u-i \ ( l - ^y-^^'(x)w(x)P u (x; I) 



W{X) 

+ a l - u (x)w'(x)V u (x; I) + a(x) l - v w(x)V ' v {x; I) 

= 7 \ f, - 1 Hxy-^wixW - v)a\x)V u 
w{x)dx l ~ u - xK v ' v nv ' v ; 

+ (t-<j , {x))V v (x]1)+<jV' v (x]1)]) 
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Comparing the rhs of Eq. proves Eq. §5§ by induction along with the re- 
cursive ODE (JHD which allows constructing systematically the complementary 
polynomials starting from Vq(x; I) = 1. For example, v — of the recursive 
ODE ® confirms Eq. ©. 

In terms of a generalized Rodrigues representation we have 
Theorem 2. The polynomials V u (x; I) satisfy the Rodrigues formulas 

d v 

V u {x; I) = w~\x)a u ~ l {x) — [w{x)a l {x)]; (10) 



V v (x; I) = w- 1 (x)^- l (x)^- JI (a l -^x)w{x)V^x- 1)) . (11) 

Proof. We prove the Rodrigues formulas for the V u (x;l) polynomials 
by integrating first the homogeneous ODE ([6]) while dropping the inhomo- 
geneous term V u+ i(x; I). This yields 

\nV u (x; I) - In c v = (—l + u + 1) \na(x) - a/ / -^c ~k [ ~~7~T> (12) 

j cr(xj J a(x) 

where c v is an integration constant and / ^fy,/ are indefinite integrals. 
Exponentiating Eq. (fT2l) we obtain 

P^(z;/) = c !/( r(x)^ +1 e" a ^^"''^^I. (13) 

Note that, if the zeros of cr(x) are real, they lie outside the fundamental 
interval (a, b) of w(x) and the hypergeometric Eq. (0Q) by definition, while x 
lies within it. So, these zeros pose no problem for the indefinite integrals. 

Now we allow for the x dependence of c„ and vary it to include the 
inhomogeneous term V u+ i(x;l). Differentiating Eq. (fl3|) and substituting 
the recursive ODE (J6]) yields 

V u+1 (x;l) = c' u (x)a- l+u+2 {x)e- ai f ^)- bl f ^) , (14) 



or 



V u+1 (x;l)a l - 2 - u (x)e ai f^ + f^) = —[ a l -"-\ x )e ai S I ^)V„(x; I)] 

= Cy{x). (15) 
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Noting that the expression in brackets on the rhs of Eq. (JT5|) differs from the 
coefficient of V u +i(x', I) on the lhs only by one unit in the exponent of a(x) 
suggests iterating the differentiation and then replacing u + 1 — ► v. This leads 
to the formula 

V u (x; I) = a- l+1+ »e- a < f^f^£- [^"V f ^ I (16) 
Integrating Pearson's ODE (j4]), 

\nw(x) = ( \ ] dx = — lncr(x) + ai f —-— + h [ (17) 

J \<j o ) J a(x) J a(x) 

and exponentiating this, gives 

i F dx | i f xdx 

W (x) = (7 6 l (18) 

Substituting this result into Eq. (fl6|) allows casting it in the general form of 
Eq. (HO]). 

When we carry out the innermost differentiation in Eq. (flOl) we obtain 
the first step (/i = 1) of the inductive proof of the generalized Rodrigues 
representation of Eq. (TTTi) . Equation (TTTi) yields trivially V v {x\ I) for [L — v, 
while for fx = v— 1 it reproduces Eq. (jSJ) and the case fx — 1 is Eq. ffTOl) . The 
inductive step from // to /i + 1 is similar to that leading to Eqs. (jSJ) and (jSJ). 

3 Generating Function 

Definition. The generating function for the polynomials V v {x\ V) is 

oo V 

V(y,x;l) = J2^'P u (x;l). (19) 

The series converges for \y\ < e for some e > and can be summed in closed 
form if the generating function is regular at the point x. 

Theorem 3. The generating function for the polynomials V v [x\ I) is 
given in closed form by 

V(y, x; I) = W( "7; (X)) f < X + V< X )) ) 1 ; (20) 
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<9 M IN w(x + wcr(x)) /a(x + ya(x))\ l M _ , . . . . , 

-7>(y,x;0 = — r~r~~^ \ r^ 21 ) V»(x + ya(x);l). (21) 



dy 11 '* w(x) V a(x) 



Proof. Equation ( 1201) follows by substituting the Rodrigues representation, 
Eq. (TTOT) in Eq. ( fT9i) which yields, with 2; = x + ycr(x), 

V(y,x;l) = f]^(x;/) = (^(xy(x))" 1 f:^^^ ; (a / (x)«;(x) 
w(x)a l (x))' £ K —^j- v ^\z)w{z)) \ z=x , (22) 



y=0 



converging for |t/cr(x)| < e for a suitable e > if x G (a, 6) is a regular 
point of the generating function, i.e. w is regular at x and x + ya(x). The 
series can be summed exactly because the expression inside the derivatives is 
independent of the summation index v and we deal with the Taylor expansion 
of the function a l (z)w(z) at the point x with increment ya{x). 

Differentiating Eq. ( fT9j) and substituting the generalized Rodrigues for- 
mula ( TTTT) in this yields Eq. ( 12TI) similarly. 

In preparation for recursion relations we translate the case /i = 1 of 
Eq. (l2~Tj) into partial differential equations (PDEs). 

Theorem 4. The generating function satisfies the partial differential equa- 
tions (PDEs) 

(1 + ya\x) + ^y 2 a"a(x)) dV{y dy X;l) = [V^x; I) + ya(x)V[(x; /)] 

• V(y,x;l), (23) 

dV{y,x;l) 



dy 



[{I - l)a'{x + ya{x)) + r(x + ya{x))}V{y, x;l-l), (24) 



1 + ya (x) + -y a a{x) 1 



2 a y ' J dx 



= V(y, x- Z)y|(l + ya'(x))V[(x; I) - l -ya"V x {x- /)}, (25) 

a{x) dV{V dx X]l) = ( 1 + y°'( x ))Mx) + (i-iW(x) + y^)(T' + (i-iW')] 

■ P(y,x;l-l)-[r(x) + (l-iy(x)]V(y,x;l). (26) 
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Proof. From Eq. (I2ip for /i = 1 in conjunction with Eq. (J5DJ) we obtain 



following from Eq. (jSJ), we verify Eq. (1231) . Using the exponent I — 1 instead 
of I of the generating function we can similarly derive Eq. ^2M . 

By differentiation of the generating function, Eq. fl22|) . with respect to 
the variable x we find Eq. (1261) . Using the exponent I instead of Z — 1 of the 
generating function in conjunction with Eq. ( 1261) leads to Eq. ( 1231) . 

4 Recursion and Other Relations 

Our next goal is to rewrite various PDEs for the generating function in terms 

of recursions for the complementary polynomials. 

Theorem 5. The polynomials V v {x\ I) satisfy the recursion relations 



a(x + ya(x)) ^ — - — = a(x)[r(x + ya(x)) 




+ (l-l)a'(x + ya(x))]V(y,x;l). (27) 
Substituting in Eq. (127]) the Taylor series-type expansions 



a(x + ya{x)) = a(x)(l + ya'(x) + -y 2 a"a(x)), 
a'(x + ya(x)) = a'(x) +ya"a(x), 

t(x + y<j(x)) = t(x) + yr'(x)a(x) (28) 



V v +x{x\ I) 



V u+1 (x;l) 



+ 



[T{x) + {l-l-v)a'{x)YP v {x;l) 

va{x)[r' + l))a"\P u ^{x- 1); 

[r(x) + l)a'(x)}V u (x; I - 1) 



(29) 



+ 




vg{x)[t' + {l-l)a"YP v -i{x;l-l) 
V x {x\ l)V v {x- l-l) + V[(x;l- l)V v -i(x] I - 1); 



(30) 
dV v {x\ V) 



dx 



(31) 
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Proof. Substituting Eq. (lT9"j) denning the generating function in Eq. (1251) 
we rewrite the PDE as Eq. (1291) . The recursion (1301) is derived similarly from 
Eq. (|24|) . The same way Eq. (1251) translates into the differential recursion 
relation (l3~Tj) . 

Corollary. Comparing the recursion (129]) with the recursive ODE ([6]) we 
establish the basic recursive ODE 

^PM = yV + {l-l-\{v- l)a")]?Vi(x; I) (32) 

with a coefficient that is independent of the variable x. 
Parameter Addition Theorem. 

V(y, x- h + l 2 )V(y, x- 0) = V{y, z; h)V{y, x; l 2 ). (33) 

E(j [^(^;^i+W- M (a;;0)-^(z;Z 1 )P^(a;;i 2 )] = 0. (34) 

Proof. The multiplicative structure of the generating function of Eq. ( 1201 ) 
involving the parameter / in the exponent implies the identity (|33l) . Substi- 
tuting Eq. (|T9|) into this identity leads to Eq. (|34|) . 

We can also separate the / dependence in the polynomials using Eq. (|28|) 
in the generating function, Eq. (12 ip . If a(x) = constant (as is the case for 
Hermite polynomials), the generating function only depends on the weight 
function, and the Taylor expansion of w (x + ya(x)) for ||/(t(x)| < 1 is equiv- 
alent to the Rodrigues formula (|T0|) . 

Corollary 1. 

<jj i n w{x + ya{x)) i 

V(y,x;l) = (l+ya(x)) 

w(x) 

- £ ! '\-£ SK (--J ww ^Sl' (35) 



— P m (x;0y(x) N - m (36) 
ml 

Proof. When a'(x) ^ 0, the Taylor expansion of the weight function in 
conjunction with a binomial expansion of the Ith power of Eq. f[2~8"j) yields 



N-l<m<N \ iV "V 
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Eq. (j35p . Using Eq. (I19p this translates into the polynomial expansion (IHoT) 
that separates the complementary polynomials Vn{x', into the simpler 
polynomials V m (x] 0) and the remaining a'(x) and I dependence. Pearson's 
ODE (jl]) guarantees the polynomial character of the V m {x\ 0) that are defined 
in Eq. (TTOj) . 

Let us also mention the following symmetry relations. 
Corollary 2. If cr(— x) = (— l) m o(x), w(— x) = (—l) n w(x) hold with integers 
m,n then Pi(-x) = (-l) l(m+1) P t (x) and 

V v {—x\V) = V v (x;l), m odd, (37) 
V v {-x;l) = {-VfV v {x-l), m even. (38) 

Proof. The parity relation for Pi(x) follows from substituting — x in the 
Rodrigues formula (J3]). The other polynomial parity relations follow from 
the identities 

V(y,—x;l) = V(y,x;l), m odd (39) 
V(—y,—x;l) = V(y,x;l), m even (40) 

which, in turn, result from substituting — x into the first formula of Theorem 
3. Expanding the generating functions according to their definition yields 
the relations (l37jl , (j55jl. 



5 Sturm— Liouville ODE 

Theorem 6. The polynomials V v [x\ I) satisfy the Sturm-Liouville differen- 
tial equation 

^ (a(x) l - u+ Mx) ^j^ ) = -\uCr(x) l -»w(x)V u (x; I), (41) 
which is equivalent to 

G{x) d2V dx^ l) + [{l ~ U)G ' {X) + T(x)] ^d^ = - KV ^ 0. ( 42 ) 
and the eigenvalues are given by 

A„ = -u[(l - ^r^y + A v = 0, 1, . . . . (43) 
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Proof. This is derived by natural induction again. The first step for v = 1 
is straightforward to verify. 

The step from v to v + 1 proceeds from the lhs of Eq. (j4"Tl) for v + 1, 
where we replace by 7^ using the recursive ODE (Q so that, after some 
elementary manipulations, we end up with 



4- ( <r(x) l -»w(x) d ' Pv+ ^ X ' ,l A = au-y ' '„■(,•) 
ax \ ax 



+ [(Z - v - la\x) + r(x))][a(x)^^ + [{I - !/)«/(*) + r(x)]^; 

+ a (x)^-[a(x) ^^ + [(/ - zz)a'(x) + r(x)] ^ i0 ; 

= a(x) l ' u - l w(x){[(l -v- l)a" + r']^ +1 (x; I) - X u V u+ i(x; I)} 

= -\ u +iPu+i(x;l), (44) 



where we have used the recursive ODE ([6]) and the ODE ( 1421) for the index 
v repeatedly. Eq. ([6]) introduces a third derivative of V v {x\ I), a term which 
shows up as the next to last term on the rhs of the first equality sign in 
Eq. (|44p . This completes the proof by induction and establishes the recursion 



Ah-i = X v -[(l-u- l)a" + r') (45) 
for the eigenvalues, whose solution is Eq. (j4"51) . 



6 Classical Polynomial Examples 

In the case of Hermite polynomials [1] , [5] , [6] , [7] , [8] a has no roots, so a(x) = 
constant = 1, without loss of generality, and a' = 0; moreover, we may 
take ai = 0, b\ = —2 so t(x) = —2x. Hence Pearson's ODE yields the 
weight function w(x) = e~ x on (— oo, oo) that is characteristic of Her- 
mite polynomials. The Rodrigues formula ([3]) then identifies the polynomials 
Pi(x) = (—l) l Hi(x) as Hermite's, while the Rodrigues formula (fTUI) for the 
complementary polynomials implies V v {x\X) = P v (x), so they are indepen- 
dent of the index I and also Hermite polynomials. The recursive ODE ([6]) 
becomes the well known differential recursion H n+ i(x) = 2xH n (x) — H' n [x). 
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The Sturm-Liouville ODE becomes the usual ODE of the Hermite polyno- 
mials. The recursion ([29]) is the basic H n+ %(x) = 2xH n (x) — 2nH n _ 1 (x). 
Eq. (132|) gives the differential recursion H' n (x) = 2nH n _i(x). The parity re- 
lation is also the well known one. The generating function is the standard 
one. Equation fl35l) reproduces the usual expansion of Hermite polynomials 
in powers of the variable x. 

For Laguerre polynomials, a(x) has one real root, so a(x) = x and 
t{x) = 1 — x without loss of generality. Pearson's ODE gives the famil- 
iar weight function w(x) = e~ x on [0,oo). Rodrigues formula ([3]) identifies 
P;(x) = l\Li{x). The Sturm-Liouville ODE (jl2|) 

x f^lll + + x) dVu ^ l) = -X v Vu(x] I), K+x = X v + 1 (46) 

allows identifying V v {x\ I) = v\L l ~ v {x) as an associated Laguerre polynomial. 
So, in the following we shift / — > I + u, as a rule. The recursive ODE ([6]) 
yields the differential recursion 

dT l (r) 

[y + l)L l - + \(x) = x^l + x)L l u (x) (47) 

which, in conjunction with 

L l - + \(x) = Li +1 (x) - L l u (x), (48) 
leads to the standard three-term recursion 

dji ( r \ 

{v + l)L l v+1 {x) = {I + v + 1 - x)L l u (x) + x^^. (49) 

The formula fflUj) of Theorem 2 is the usual Rodrigues formula for associated 
Laguerre polynomials, while the generalized Rodrigues formula (fTTj) 



LUx) = ^e x x~ 1 ^- ( x l+u -»e- x L l + u -»(x)) (50) 



is not part of the standard lore. 

The generating function ( J20l) for this case becomes 



oo 



L(y, x- 1) = £ rt-'Cs) = e~ xy (l + y) 1 (51) 



u=0 
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and is simpler than the usual one for associated Laguerre polynomials, which 
is the reason why our method is more elementary and faster than the standard 
approaches. The recursion (129]) becomes 

(u + l)L l - + \(x) = (I - x)Ll(x) - xL l +_\(x), (52) 

while the recursion fl30l) becomes 

{y + l)L l v+l (x) = (l + u + l- x)L l u (x) - xL l +\{x), (53) 

and Eq. (1311) translates into 

+ = _ L m ( x) - Lltl(x), (54) 

a sum of the known recursion = —L^^x) which is the basic recursive 

ODE (I3"i2"|) . Equation (1551) gives the standard expansion 

The simplest addition theorem originates from the elegant identity 

L(y, xi, ni)L(y, x 2 ; n 2 ) = L(y, x x + x 2 ; n x + n 2 ) (56) 
which translates into the polynomial addition theorem 

V v {xi + x 2 ; m + n 2 ) = ( i ) 'Pv-k{xx\ ni)V k {x 2 ; n 2 ) (57) 

fc=0 Vfc/ 

and 



- ni+n2 



[x l + x 2 ) = Y J LT-\x 1 )L n u l\\x 2 ) (58) 



k=0 



for associated Laguerre polynomials which is not listed in the standard ref. j8] 
or elsewhere. 

In the case of Jacobi polynomials, a(x) has two real roots at ±1, without 
loss of generality; so 

a(x) =(l-x)(l + x), t(x) = b-a-(2 + a + b)x, (59) 
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in a notation that will allow us to use the standard parameters. 
Pearson's ODE (0} leads to 

w(x) = (l-x) a (l + x) b , (60) 

and Rodrigues formula ([3]) and ({TO]) identify the polynomials 

P^x) = 2\-l) l l\Pi a ' b \x), V u (x; I) = (-2Yv\Pi a+l - v > b+l - u \x). (61) 

Thus, we shift I — > I + v in translating our general results to Jacobi polyno- 
mials, as a rule. We may also set / = because this index merely shifts the 
parameters a, b. 

The recursive ODE ([6]) translates into 

-2(iy+l)Pj / a +1 hb ' 1 \x) = \b-a-(a + b)x]P^ b \x) 

f ip(a,b)( \ 

+ (i- x y^ [x) . (62) 

ax 

The Sturm-Liouville ODE (142]) reproduces the usual ODE of Jacobi poly- 
nomials. The generating function, Eq. ( 1201 . 

[1 - x - y(l - x 2 )] a [l +x + - x 2 )] 6 f 1 - (x + y(l - x 2 ) 2 ~\ l 
} ~ (l-x) a (l+x) 6 I l^x~ 2 J 

= [l-y(l + x)} a [l + y(l-x)] b [l-2xy-y 2 (l-x 2 )} 1 (63) 

is much simpler than the standard one [I], especially when we set / = 0, allow- 
ing for the transparent derivation of many recursion relations. For example, 
Eq. (1301) becomes 

-4(z/+l)P^ 1 '^ 1) (x) = 2[b-a-x(a + b)]Pj ) a ' b \x) (64) 

- (l-x 2 )[v + l + a + b]P^ h+1 \x), 

Eq. (130]) translates into 

- 4(i/ + 1) Pfrx'^ix) = 2[b-a-x(a + b + 2u)]Pj ) a - 1 ' b - 1 \x) 

+ (1 -x 2 )[a + b + 2u]P^ b \x), (65) 

and Eq. (|3~T1) takes the form 

(x2 _ 1} ^a" b+2) (x) + ix dP^ b+i \x) + 4 _dp^)(x) 



doc doc doc 

2[a + b + 2v]Pl a+ 1 1,b+1 \x) + [b-a + x(a + b + 2v)\ 

i>+ 2 ' 6+2) (x). (66) 
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Equation (155|) gives 

TV , x 



P^ b \x) = (-2) N N\(1 -x)-(l + z)-> ( ) 

A product formula for Jacobi polynomials is obtained from an addition 
theorem in the variable y for our generating function for I = (where we 
display the upper parameters now for clarity) 

V^ b \ yi , x- 0)V {a > b \y 2 , x- 0) = [(1 + yi (l -x)(l + y 2 (l - x)] b 
■ [(l- ?/l (l + x)(l-2(l + x)] a 

= [1 + (m+y2)(1 -, )r{1 + _^gL_ r 



V^ b \ yi + y 2 ,x;0) Y, 

j,k=0 



a\ fb\ ( 2/i|/2(l - x 



,2 



kJ \jj \l + (yi + 2/ 2 )(l - x) 



\1 - (yi + y 2 )(l + x)J jj^o\kJ\jJ 
. V (a-k, b - 3){yi + y2) x . )^+*yj+ fc (i - x) 2 ^(l + x) 2k . (68) 

Expanding into Jacobi polynomials according to Eq. (1191) . comparing like 
powers of yiy 2 , converting to Jacobi polynomials and shifting a — * 
a + Pi , b — > b + v\ yields the product formula 

p( a,fc ) (x)P( aJrV1 ~ V2,bJrl/1 ~ V2 \x) = Y 2~ 2v ( Vl + V2 ~ ^ 



0<v<(^i+y 2 )/2 

2fc n _ ^ 2 O-*0 



1/1 — 1/ 



fc=0 

• 4+i 1/ -2 !y r fe ' b+ ^ 2+fe) w- (69) 

7 Conclusions 

We have used a natural way of working with the Rodrigues formula of a given 
set of orthogonal polynomials which leads to a set of closely related comple- 
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mentary polynomials that obey their own Rodrigues formulas, always have a 
generating function that can be summed in closed form leading to a transpar- 
ent derivation of numerous recursion relations and addition theorems. These 
complementary polynomials satisfy a homogeneous second-order differential 
equation similar to that of the original polynomials. 

Our method generates all the basics of the Hermite polynomials. It gen- 
erates the associated Laguerre polynomials and many of their known prop- 
erties and new ones from the Laguerre polynomials in an elementary way. It 
also simplifies the derivations of various results for Jacobi polynomials. 

Our method is not restricted to the classical polynomials; when it is ap- 
plied to the polynomials that are part of the wave functions of the Schrodinger 
equation with the Rosen-Morse and the Scarf potentials, it links these poly- 
nomials to the Romanovski polynomials which will be shown elsewhere. 
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